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The observation of light nuclei at ALICE and the X(3872) conundrum 
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The new data reported by ALICE on the production of light nuclei with p±^ < 10 GeV in Pb-Pb 
collisions at ^snn = 2.76 TeV are nsed to compute an order-of-magnitude estimate of the expected 
production cross sections of light nuclei in proton-proton collisions at high transverse momenta. We 
compare the hypertriton, helium-3 and deuteron production cross sections to that of A(3872), mea¬ 
sured in prompt pp collisions by CMS. The results we find suggest a different production mechanism 
for the A(3872), making questionable any loosely bound molecule interpretation. 
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As first discussed in [1], one expects a suppression of 
loosely bound hadron molecules in high energy pp{p) col¬ 
lisions. Small relative momenta in the center of mass 
of such molecular hadrons, needed to preserve a state 
with few keVs’ binding energies, are in fact hard to be 
obtained in hadron collisions at high energy and p±. 

Despite this, the A(3872), one of the most studied 
loosely bound hadron molecule candidates [2] , is strongly 
produced at the LHC - see e.g. [3]- This might simply 
be the indication that the X hadron molecule interpreta¬ 
tion is not correct (for the alternative tetraquark model, 
see |1HS])- 

Assuming that final state interaction mechanisms are 
at work - whose description requires several model- 
dependent assumptions 13 [g - it has been proposed that 
the relative kinetic energy might be reduced in the cen¬ 
ter of mass of the hadron pair constituting the X, in 
such a way to match a shallow discrete level of some 
inter-hadron potential. A hadron molecule would then 
be formed, with a precise relation between binding en¬ 
ergy and strong coupling to its constituent hadrons [S]. 
Since the mass and branching ratios of the X have not 
been measured with the required precision yet, it is still 
unclear if this relation is fulfilled. 

Final state interactions should also favor the prompt 
formation of bona fide light nuclei in high energy hadronic 
collisions. It would therefore be of great interest to mea¬ 
sure the pp (anti) deuteron production cross section in the 
same p± region where the X has been observed [10]. 

Unfortunately, (anti)deuteron production in pp colli¬ 
sions at p± values as high as « 15 GeV (where the X is 
clearly seen at CMS 0 ) has not been measured yet. 

However, very recently the ALIGE collaboration re¬ 
ported results on the production of deuteron, helium-3 
(^He) and hypertriton (^H) light nuclei in relatively high 
p± bins in Pb-Pb collisions, at yfsVN = 2.76 TeV [TTlITg . 
This is potentially a very exciting result for the reasons 
described above. 

We would like to draw the attention on these data and 


propose a way to exploit them to provide an order-of- 
magnitude estimate of light nuclei production in pp col¬ 
lisions, to compare with the X data. 

As a first approximation one can assume that there are 
no medium effects enhancing or suppressing the produc¬ 
tion of light nuclei in Pb-Pb collisions. This is equivalent 
to state that each nucleus-nucleus collision is just an in¬ 
dependent product of A/ou proton-proton collisions, with 
fVcoii computed in a Glauber Monte Carlo calculation as 
a function of the centrality class. We use the results 
from m, which are compatible at Itr level with the AL¬ 
ICE ones [13], and never more different than 3%. To 
compare with ^/s = 7 TeV data, we rescale our estimated 
cross sections by a factor cr“®'(7 TeV)/cr“®'(2.76 TeV) = 
1 . 1 . 

Consider for example the production of hypertriton 
observed by ALICE in Pb-Pb collisions Neglecting 
medium effects, the pp cross section can be estimated 
with 

dcT (jH) / 

dp± / 

/ pp 

Ay 1 / d^jV(3He7r) \ 

6(3He7r) ^ Cpp V dpj_dy J 

Ay / d^iV(3He^) \ 

;B(3He7r) ^ Vevt V dp±dy ) 
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ALICE analyzes ^Heyr pairs, thus we need to divide 
by the branching ratio for the —>■ ^HeTr decay - 

5(^He7r) 25% [T^ - in order to deduce the number 


^ In the following, the average of hypertriton and anti-hypertriton 
data is understood. 
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FIG. 1: Comparison between the prompt production cross section in pp collisions of X(3872) (red), deuteron (green), ^He 
(orange), and hypertriton (blue). The X data from CMS [3] are rescaled by the branching ratio B{X —>■ J/‘4nni). Deuteron 
data in pp collisions are taken from ALICE m- The ®He and hypertriton data measured by ALICE in Pb-Pb collisions [III 112) 
have been rescaled to pp using a Glauber model, as explained in the text. The dashed green line is the exponential fit to the 
deuteron data points in the pa_ £ [1.7, 3.0] GeV region, whereas the dotted orange one is the fit to the ®He data points. The 
solid and dot-dashed blue lines represent the fits to hypertriton data with Raa = 1 (no medium effects) and an hypothetical 
constant value of Raa = 5. The hypertriton data points are horizontally shifted at the bin centres of gravity - being defined 
as the point at which the value of the fitted function equals the mean value of the function in the bin. (Left Panel) The 
hypertriton data are fitted with an exponential curve, and the light blue band is the 68% C.L. for the extrapolated Raa ~ 1 
curve. ®He data in the p± € [4.45, 6.95] GeV region are also fitted with an exponential curve. (Right Panel) The hypertriton 
and ®He data are fitted with blast-wave functions m whose parameters are locked to the ®He ones obtained in m- 


of parent hypertritons. We stress that the experimen¬ 
tal data in m are indeed normalized to Nevt = -^p'b.pbJ 
i.e. the total number of inelastic Pb-Pb collisions anal¬ 
ysed (about 20 X 10® events in the 0-10% centrality bin). 
We use cr™®' = 73 mb, as measured in -/s = 7 TeV colli¬ 
sions [TB], and Ay = 2.4 to compare with the CMS analy¬ 
sis [3]. In this centrality class, we use = 1518 [13]. 


Similarly, we can estimate the ®He distribution in pp 
collisions from the ALICE Pb-Pb data in the 0-20% cen¬ 
trality class [13], using = 1226 [TS]. We re¬ 

mark that the selection of these events rejects any ®He 
not produced in the primary vertex, i.e. the hypertri¬ 
ton decay products. Since the ®He data points with 
p± < 4.4 GeV show a deviation from the exponential 
behavior, likely due to the expansion of the medium, we 
perform an exponential fit to the points in the region 
p± G [4.45, 6.95] GeV only. Alternatively, we fit hypertri¬ 


ton and ^He data with the blast-wave model which de¬ 
scribes particle production properties by assuming ther¬ 
mal emission from an expanding source m- This model 
is expected to reproduce correctly the low and medium 
p± regions in Pb-Pb collisions. Since we are rescaling 
Pb-Pb data to pp by a constant factor, the same shape 
holds in our estimated pp data, and gives a guess on the 
asymptotic exponential behavior. The results are shown 
in Fig. [2 

Our rescaling to pp collisions does not take into ac¬ 
count neither medium effects, nor the fact that the co¬ 
alescence/recombination mechanism can be enhanced in 
Pb-Pb collisions m- In fact, such phenomena are known 
to favor the production of many-body hadrons with re- 


^ The blast-wave function is 


dN 

dpx_ 


ocpj_ 


l\drm^lo[P^^)K, 


m± cosh p 
^kin 


where is the transverse mass, R is the radius of the fireball, Iq 
and are the Bessel functions, p = tanh“^ ^ (^+2)(^) ^ 

and (/3) the averaged speed of the particles in the medium. 
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FIG. 2: Comparison between the nuclear modification fac¬ 
tor Rcp for deuteron (red) and for generic charged tracks 
(blue) |23| in central (resp. 0-10% and 0-5%) versus periph¬ 
eral (60-80%) Pb-Pb collisions. We evaluate Rcp either tak¬ 
ing the bin-by-bin ratio of ALICE data [12] (full disks), or 
using the values given by the blast-wave fit for the periph¬ 
eral Pb-Pb data, and extrapolated up to p±_ ~ 4.8 GeV (open 
disks). The dashed line corresponds to no medium effects, 
Rcp = 1. It is worth noticing that Rcp{d) gets enhanced at 
p± > 2.5 GeV. 


sped to what expected in vacuum. Medium effects are 
discussed later. 

For the deuteron we use ALICE pp data [12] to esti¬ 
mate 


f da{d) \ 1 d^Nid) \ 

Wpu j,, ^ ^^[N^rdp±dy) 


( 2 ) 


being the number of pp inelastic collisions col¬ 
lected. We perform the fit to the points in the region 
p± S [1.7,3.0] GeV, which shows a good exponential be¬ 
haviour. 

The CMS analysis of X production provides the 
differential cross section times the branching fraction 
B {X{3872) —>■ J/i/; 7 r+ 7 r“). The latter has not been mea¬ 
sured yet, and the lower limit reported in the PDG is 
B > 2.6% [Hj. An estimate for the upper limit has 
been reported, B < 6 . 6 % at 90% C.L. [23; we use in¬ 
stead the more conservative value B = 8 . 1 + 3 ;®% [H]. The 
comparison in Fig. shows that, according to the most 
conservative exponential fit in the left panel, the extrapo¬ 
lated hypertriton production cross section in pp collisions 
would fall short by about 2 4- 3 orders of magnitude with 
respect to the X production, and much more according 
to the blast-wave fit in the right panel. The drop of the 
deuteron cross section, which is directly measured in pp 
collisions, appears definitely faster. 

As we mentioned already, the main problem for the 
production of loosely bound molecular states in proton- 
proton collisions is the difficulty in producing the con¬ 
stituents close enough in phase space. However, it is 
well known that the interaction of elementary partons 


with the collective hot dense medium causes relevant en¬ 
ergy loss of the partons themselves. This effect is usually 
quantified by the nuclear modification factor pTVBS] 


Raa 


d^N 




^ A^evt dp±dy j 

d‘^N \ 


R7co\\ ( 


Wvt dp±dy J 
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which compares the particle yield in Pb-Pb collisions 
with that in pp. It then follows that the method used 
to obtain Eq. 0 corresponds to assume Raa = 1- 

While for ordinary hadrons medium effects generally 
lead to a suppression of the particle yield - i. e. Raa < 1 
- conversely they can favor the production of hadronic 
molecules. The role of the medium would be, in fact, that 
of decreasing the relative momenta of the components 
with respect to the zero temperature case due to the well- 
known jet quenching effect [21 HTj. This would favor 
their coalescence into the final bound state by reducing 
their relative momenta directly at parton level. 

The coalescence model is based on the sudden approx¬ 
imation ® and is implemented by calculating the overlap 
of the density matrix of the constituents with the Wigner 
function of the final composite particle. In particular, it 
has the important property of taking into account the 
inner structure of the considered hadron. If one only re¬ 
quires vicinity in momentum space, the p± distribution 
of a composite state with N constituents coming out of 
a hot QCD medium is roughly given by 


dN^ 

dp± 


N 


ip±) - n 


dN, 

dp± 


(PJN), 


(4) 


where W is the number of final bound states and W is 
the number of produced constituents. This would also 
explain why in Fig. the cross section for the ®He and 
hypertriton are several orders of magnitude smaller than 
the deuteron one: one additional p or A, close enough in 
phase space, must be produced. 

It has already been shown that coalescence effects in 
Pb-Pb collisions can have relevant consequences on the 
production of multi-quark states. In particular, molecu¬ 
lar states with small binding energy are expected to be 
enhanced, i.e. Raa > 1 [28] . 

Unfortunately there is no measurement of Raa for the 
deuteron as a function of p± . However, there is another 
nuclear modification factor which is often used, 

( 1 d^N 

\ Wvt dp±dy 
Hcp = -™™ 

( 1 ' 

I Wevt dp^dy I 


) U-iU70 j 

Pb-Pb / 


ArO-10% 


DU-5U70 f 
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i.e. the assumption that the binding of the constituents happens 
on small time scales and therefore their wave function remains 
unchanged during the transition to the bound state. 
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This quantity is a comparison between the most cen¬ 
tral and the most peripheral Pb-Pb collisions and there¬ 
fore provides another valid indicator of the strength of 
medium effects (which should be absent in the less dense, 
most peripheral events). The fact that Raa and Rcp 
measurements for hadron species are strongly correlated 
to each other is shown experimentally by a thorough 
data analysis reported by ATLAS |23], up to very high 
p_L ~ 100 GeV. 

Using the ALICE data presented in m we can com¬ 
pute Rcp for deuteron as a function of p± and compare 
it with that for generic charged tracks, as reported in [23] 
- see Fig. We use = 27.5 [13. As one imme¬ 

diately notices, the difference from ordinary hadrons is 
striking. The presence of the QCD medium is extremely 
effective at enhancing the production of deuteron for the 
reasons explained before. In fact, Rcp for this hadronic 
molecule becomes larger than unity for p± > 2.5 GeV, 
in particular we have Rcp = 1.7 at the last point with 
p± = 3.1 GeV. Using the blast-wave fitting function for 
the peripheral data taken from [IB, we also extrapolate 
up to the end point of the central data, confirming the 
growth of Rcp with p±. 

We expect a similar behavior in Raa^ in particular a 
value larger than 1 for p± large enough. 

To get an independent rough estimate for Raa, we 
assume the deuteron production cross section in pp col¬ 
lisions to scale with ^/s like the inelastic cross section, 
and compare the ALIGE data in central Pb-Pb colli¬ 
sions at = 2.76 TeV, with the ones in pp collisions 

at y/s = 7 TeV [12]. Indeed, we find that Raa exceeds 1 
at p± = 2.1 GeV, and reaches 5 at = 4.3 GeV. This 
gives strength to our expectation for Raa > 1. To dis¬ 
play the size of this effect, we plot also the hypertriton 
curves for Raa = 5 in Fig. [T] 

One naturally expects for a similar enhancement to 
be even more relevant for 3-body nuclei like ^He and 
hypertriton. Its role would be to further decrease the 


extrapolated cross section in prompt pp collisions. As 
we already said, indeed, a value of Raa > 1 applied to 
Pb-Pb data implies a pp cross section even smaller than 
predicted by the Glauber model. Even though qualitative 
conclusions can already be drawn, a quantitative analysis 
substantiated by data at higher p± is necessary for a 
definitive comparison with the X case. 

Even assuming that only a hot pion gas is excited in 
Pb-Pb collisions, there would likely be a large number of 
final state interactions with pions catalizing the forma¬ 
tion of a loosely bound hypertriton along the lines dis¬ 
cussed in [ittniiis]. In any case, such an environment 
is present in the Hadron Resonance Gas corona formed 
when the outer shell of the QGD medium cools down [30] . 

In summary, the extrapolation of deuteron and ^He 
data in pp collisions shown in Fig. [^ suggests that loosely 
bound molecules are hardly produced at high p±. The 
extrapolated curve of hypertriton data from Pb-Pb col¬ 
lisions might lead to milder conclusions although we ex¬ 
pect it should be significantly suppressed when medium 
effects are properly subtracted. Such effects are indeed 
already sizeable for the deuteron as shown in Fig. [^ and 
probably even more relevant for 3-body nuclei. 

We are aware that for an unbiased and definitive com¬ 
parison with X production at p± as high as 15 GeV, 
deuteron (or hypertriton) should be searched in pp colli¬ 
sions rather than in Pb-Pb to avoid the complicacies of 
subtracting medium effects. These analyses can be per¬ 
formed by ALIGE and LHGb during Run 11. One of the 
purposes of this letter is to further motivate the required 
experimental work. 
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